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ABSTRACT 


Context. The chemical composition of planets is an important constraint for planet formation and subsequent differ¬ 
entiation. While theoretical studies try to derive the compositions of planets from planet formation models in order 
to link the composition and formation process of planets, other studies assume that the elemental ratios in the formed 
planet and in the host star are the same. 

Aims. Using a chemical model combined with a planet formation model, we aim to link the composition of stars with 
solar mass and luminosity with the composition of the hosted planets. For this purpose, we study the three most 
important elemental ratios that control the internal structure of a planet: Fe/Si, Mg/Si, and C/O. 

Methods. A set of 18 different observed stellar compositions was used to cover a wide range of these elemental ratios. 
The Gibbs energy minimization assumption was used to derive the composition of planets, taking stellar abundances 
as proxies for nebular abundances, and to generate planets in a self-consistent planet formation model. We computed 
the elemental ratios Fe/Si, Mg/Si and C/O in three types of planets (rocky, icy, and giant planets) formed in different 
protoplanetary discs, and compared them to stellar abundances. 

Results. We show that the elemental ratios Mg/Si and Fe/Si in planets are essentially identical to those in the star. 
Some deviations are shown for planets that formed in specific regions of the disc, but the relationship remains valid 
within the ranges encompassed in our study. The C/O ratio shows only a very weak dependence on the stellar value. 


1. Introduction 


The composition of planets is important for an understand¬ 
ing of the planet formation process, and the evolution and 
structure of planets. While the composition of the atmo¬ 
sphere of planets can be observed through the absorption 
of the star light, for example, the composition of the solids 
remains difficult to estimate. There are two main hypothe¬ 
ses used to estimate the rocky composition of exoplanets. 
In the first hypothesis, it is assumed that the solids formed 
in e quilibrium from the material of the solar nebula (Else r 

et al. l2012l : lMadhusudhan et al.ll2012l : lThiabaud et al.ll2014T >, 

so that equilibrium chemistry can be used as a constraint 
for modelling. The second hypothesis assumes that the el- 
emental ratios Fe/Si and Mg/Si have s tellar values (e.g., 
IValencia et al.l 12010 : 1 Wang et alJ l2013f ) and these values 
are then used to model the planetary structure. This last 
assumpti on well reproduces the fo rmation and composition 
of Mars dKhan fc (IonriolIvll2()()8l ) and Earth (I.Tavov et al.l 
120101. PREM) but not that of present-day Mercury, which 
has elemental rati os that are not so lar, possibly due to its 
collisional history dBenz et al.ll2007i l or peculiar formation 
rtbewidll 9721. Il 988IT 


Elemental ratios are important because they govern the 
distribution and formation of chemical species in the pro¬ 
toplanetary disc. Mg/Si governs the distribution of sili¬ 
cates. If the Mg/Si mass ratio is lower than 0.86, Mg 
is incorporated in orthopyroxene (MgSi0 3 ) and the re¬ 


maining Si enters other silicate minerals such as feldspars 
(CaAl 2 Si 2 0 8 , NaAlSi 3 0 8 ) or olivine (Mg 2 Si0 4 ). For Mg/Si 
ratios from 0.86 to 1.73, Mg is distributed between pyroxene 
and olivine. For Mg/Si higher than 1.73, Si is incorporated 
into olivine, and the remaining Mg enters other minerals, 
mostly oxides. 

The amount of carbides and silicates forme d in planets 
is controlled by the C/O ratio (lLarimerHl975 i). For mass 
ratios higher than 0.6, Si is expected to combine with C 
and form carbides (SiC), because the system is C-rich. If 
the C/O ratio is lower than 0.6, Si combines with O to form 
silicates, which are basic building blocks of rock-fornring 
minerals,for instance sioQ However, as a result of the 
chemical properties of C (used to form carbohydrates, 
CO, C0 2 , etc.), the gaseous C/O ratio in planets can 
deviate from the stellar value depending on different 
parameters (temperature, pressure, oxidation state, etc.) 
and processes during planet formation, including the 
initial location of formation of the planetary embryos, the 
migration path of the planet, a nd the evo lut ion o f the gas 
phas e of a protoplanetarv disc (lOberg et al.ll2011l: Ali-Di b 
et al. l2014l : lMadhusudhan et al.ll2014l : lThiabaud et al.ll2015l f. 

In the present study, we show that assuming the ele¬ 
mental Fe/Si and Mg/Si ratios in planets to be of stellar 

1 Note that the Mg/Si and C/O ranges given here are valid for 
pressures of ca. 1 x 1CU 4 bars. The pressures in protoplanetary 
discs are lower at most of the distances to the star. 
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value agrees with the formation of material by condensa¬ 
tion from the stellar nebula. To follow this assumption, 
we quantified the abundance of refractory compounds in¬ 
corporated into planets for different compositions of stellar 
nebulae using equilibrium chemistry. From these computed 
compositions, the elemental ratios Fe/Si, Mg/Si, and C/O 
in solids (minerals and ices) were derived for three types 
of planets: rocky, gas giant, and icy planets. These ratios 
were then compared to the ratios in the host star given in 
the literature, to determine whether it is valid to assume 
that ratios are equal in stars and planets. The goal of this 
study is not to predict the abundances of elements in the 
planets, but to provide relationships between the star and 
the planetary composition for stars of solar mass and lumi¬ 
nosity with different compositions. In Section 2, we recall 
the model used to determine the chemical composition of 
planets. Section 3 is dedicated to the sample of stars from 
which the composition of the stellar nebula is derived. Re¬ 
sults and discussions are presented in Section 4. Section 5 
summarizes the conclusions. 


2. Planet formation and composition model 

The planet formatioiynodeRrsedinthis_stud^sRiat - of,Al- 
ibert et al. (120131 ) (hereafter A13) and is based on the core 
accretion model. Planets are formed in a protoplanetary 
disc by accretion of solid planetesimals until they eventu¬ 
ally become massive enough to gravitationally bind some 
of the nebular gas, surrounding themselves with a tenuous 
atmosphere. 

2.1. Disc model 

The model of a proto pla netar y disc is structured around 
three modules (lAlibert et al.l l2005 il. which calculate the 
structuremidThe^vohitioii^amondrradiatedJ^e^Alibert 
et al . 2013 . an d references therein) and a fully irradiated 
( [Fouchet et al.ll2012l : iThiabaud et al.lf2014l : IMarboeuf et al.1 
l2014bT I 2 l protonlanetarv disc. The disc parameters are com¬ 
puted using a 1+1D model, which first resolves the vertical 
structure of the disc for each distance to the star to derive 
a vertically averaged viscosity. The calculation includes the 
hydrostatic equilibrium, the energy conservation, and the 
diffusion equation for radiative flux. The radial evolution 
is then computed as a function of this vertically averaged 
viscosity, photo-evaporation, and mass accretion rate by 
solving the diffusion equation. 

The initial gas surface density (in g.cm -2 ) is computed 
following 


of solid discs through a dust-to-gas ratio randomly chosen 
to lie between 0.008 and 0.1 from a list of ~1000 CORALIE 
targets (see A13 for more details) to produce 500 different 
discs, whose characteristics are presented in Table lA~2l 

The parametric surface density profile (Eq. [[]) is only 
used as an initial condition (t=0). This initial profile 
evolves with time, and the surface density profile at t!= 
0 is recomputed at each time step, following 

/f = - r j- r ( rl/2 |:( l/Srl/2 )) + Qacc + Qph, (2) 

where Q acc is an additional term related to the accretion 
by planets and Qph an additional term related to photo¬ 
evaporation. 


2.2. Planet formation 


The simulations of planet formation start with a protoplan¬ 
etary disc, composed of gas with planetesimals, and ten 
planetary embryos of lunar mass. These embryos are uni¬ 
formly distributed (in log) between 0.04 and 30 AU. The 
processes of planetesimal formation are not the topic of 
this contribution, and the model describes processes after 
their formation. Planetesimals are set to be large enough 
(~ 1km) so that their drift can be neglected. The model 
includes several processes: 


— accretion of planete simals from which grows the solid 
core of the planet (lAlibert et al.l 120051 : iFortier et al.l 

[Mill- 


accretion of gas. If the core of the planet reaches the 
critical core mass, the gas is accreted in a r unaway fash- 
ion, leading to the formation of gas giants (lAlibert et al.l 
1200511 . 

planet-disc interactions, whic h contributes to mig ration 
(type I or type III of planets (lAlibert et all 12005 : Mor- 
dasini et al. l2009al lblf. 

planet-planetesimal interactio ns, by which the pla netes- 
imal is excited by the planet (IFortier et al.ll2013l l. 
gravitational planet-planet interactions, making it pos¬ 
sible to catch a pair of planets (or more) in mean motion 
resonance, as well as pushing each other. 


The calculations are stopped once the gas disc has dissi¬ 
pated. The cumulative distribution of disc lifetimes is as¬ 
sumed t o decay expon etially with a characteristic time of 
2.5 Ma (IFortier et al.l 120131 A13). The photoevaporation 
rate is then adjusted so that the protoplanetary disc mass 
reaches 10 - 5 Mq at the selected disc lifetime, when calcu¬ 
lations are stopped. The longest lifetime of a disc is 10 Ma 
(see A13 for more details). 


E = £oQ^ 7 e^) 2 " 7 , (1) 

where ao is equal to 5.2 AU, a core is the characteristic 
scaling radius (in AU), 7 is the power index, r the distance 
to the star (in AU), and £ 0 =(2-7) A [ d -fy c 7 (in g.cm -2 ). 

ZTT Ct C ore 

The values for a core , 7 , and So are v aried for every disc an d 
are derived from the observations of I Andrews et al.l (|2010l ). 
The mass of the gas disc is inferred from the observations 

2 The irradiation is included by modifying the temperature 
boundarvcondition_at i ^he i disc ii surface Ji JolIowin^Hueso & Guil- 
lot (120050. 


2.3. Refractory composition 

The chemicaljnodeMs_simila]^Jothatj i ised^i^Thiabaud 
et al. (|20l4 120151 ): IMarboeuf et all (i2014al lbl). However, in 
contrast to these studies, the assumption that all discs have 
the same pressure profile (equal to the average pressure pro¬ 
file of all discs) for computing the chemical equilibrium is 
not usec0. The presence and abundances of refractory min¬ 
erals in planetesimals is ruled by the assumption of equilib¬ 
rium condensation from the primordial stellar nebula. The 

3 The implications of this assumption have been discussed in 
iThiabaud et alJ (I2014T) . 
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software HSC Chemistry (v.7.1) computes this equilibrium 
using a Gibbs energy minimization routine; the pressure 
and temperature profiles are provided by the planet forma¬ 
tion model for each disc. The composition of species in the 
disc is computed for the initial profiles of E, T, and P and 
is assumed not to change. As in iThiabaud et all (l2014f) . the 
refractory elements have condensed by the time the system 
has cooled below 200 K; the amount of refractory species 
at a lower temperature is thus given by their abundance at 
200 K. 

Thiabaud et al. (120141 ) showed that a model in which 
refractory organic compounds are not taken into account 
better reproduces the observation of the solar system. Con¬ 
sequently, the results shown in the present study do not 
consider the formation of such species, even though obser¬ 
vations on single objects (such as Haley’s comet) suggest a 
high carbon inventory that is due to these species, which 
could be the result of surface effects. The inclusion and 
discussions of refracto ry organic material can be found in 
IThiabaud et al.l (I2014T ). 

2.4. Volatile composition 
2.4.1. Volatile components in ices 

The chemical composition of eight volatile components 
(CO, C0 2 , N 2 , NH 3 , H 2 0, CH 3 OH, CH 4 , H 2 S) in plan¬ 
ets is determined by condensation and gas-trapping in 
water ice (e.g. clathrates) on the surface of refractory 
grains. The thermodynamic conditions of these processes 
are determined through the partial pressures of the volatile 
molecules in the gas phase of the disc and by the tempera¬ 
ture and surface density of the disc. If the partial pressure 
of a molecule is higher than the equilibrium pressure of 
condensation, then the species condenses and its amount 
at the dis tance r to the star is det ermined by the following 
equation dMousis fe Gautier! l2004f ): 


m x {r) 


Y x Y(T x ,P x ,r ) 

Yh 2 o Y(T H2 o,Ph 2 o,i") ’ 


(3) 


where is the mass ratio of the molecules x rela¬ 
tive to H 2 and Yh 2 o the mass ratio of H 2 0 relative to 
H 2 in the disc; Y,(T x ,P x ,r) and X (T H2 o, Pn 2 o,r) are 
the surface densities of the disc at a distance r to the 
star, when the molecule x (y^H 2 0) is condensed or trapped. 


2.4.2. Volatile components in the gas phase 

The volatile molecules considered in the solid phase are also 
considered in the gas phase of the disc. The computation 
of their abundances is similar to that of the ices and is 
calculated following: 


N x = Y x .N Sa , (4) 

where N x is the abundance in the gas of species X, Y x 
is the ISM abundance of species X relative to H 2 , and Nh 2 
is the abundance of I7 2 . 

In contrast to the solid phase, in which the planetesi- 
mals are assumed to be large enough not to drift, the gas 
composition evolves with time follo wing the provisio n given 
b y the theory of accretion discs by iLvnden-Bell fc Pringlel 
(I1974T ). The radial velocity of the gas u ac jv (in cm.s -1 ) due 
to mass and angular momentum conservation in the pres¬ 
ence of a viscosity v is calculated using 

= (5) 

The calculation of these abundan c es in the gas phase is 
the same as that of IThiabaud et all (2015j). In the latter, 
the importance of planet migration and of the gas phase 
evolution was demonstrated, so that these effects are also 
included in this study. The authors also showed that it is 
possible to enhance the C/O ratio of the gas disc to up to 
three times the stellar value, depending on the time and lo¬ 
cation, which makes the formation of specific objects, such 
as the enstatite chondrite parent bodies possible (see their 
Fig. 8). 

3. Sample 

To quantify the relationship between elemental ratios in 
stars of solar mass and luminosity and planets, a sample 
set of observed compositions from 18 stars was selected for 
this study. This sample set encompasses a wide range of 
Fe/Si (from 1.34 to 2.81), Mg/Si (from 0.66 to 1.85), and 
C/O (from 0.12 to 0.75) mass ratios. The central elemental 
abundances in these stars and their central elemental ratios 
are given in Table IA.1I and Table [Q 

For each star of the sample set, the composition of the 
protoplanetary disc was computed as discussed in the pre¬ 
vious sections, taking the stellar abundances as proxies for 
nebular abundances. 


In previous studies by iMarboeuf et al.l (l2014al fbj) , two 
models for which the initial CO/H 2 0 molar ratio is 0.2 or 
1 were computed, along with the presence of clathrates. 
The formation of clathrates in protoplanetary discs is un- 
certaiijf], and the change of the initial CO/H 2 0 does not sig¬ 
nificantly influence the abundances of refractory elements, 
including the Fe/Si and Mg/Si ratios. Moreover, the C/O 
ratio is not sign i ficant ly changed in planets, as shown in 
IThiabaud et al.l (120151 ). if the initial CO/H 2 0 is equal to 
0.2 or 1. For these reasons, clathrates are not included, 
and we only computed the model with an initial CO/H 2 0 
= 0 . 2 . 


4 The kinetics for clathrate formation are poorly constrained 
at low temperatures a nd pressures, and the kinetics feasibility 
remains uncertain (see IMarboeuf et all (12015!) for more details) 


4. Results and discussions 

The values for elemental ratios of planets were computed by 
distinguishing three types of planets: rocky, icy, and giant 
planets. Each type of planets forms a population, based on 
its mass and position. Rocky planets were defined to be 
planets whose mass is lower then 10 M® and whose semi¬ 
major axis is lower than 2 AU. Icy planets were defined as 
low mass planets below 10 M® located beyond 2 AU, and 
giant planets were defined as massive planets with masses 
greater than 10 M®@. The mode valuc|j of the composition 

5 In this definition, Neptune- and Uranus-type planets are not 
considered as icy planets but as giant planets. 

6 The mode value of a sample is the value that appears most 
often in a set of data. 
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Table 1: Elemental ratios of stars in the sample set, in 
mass and moles. 


Star 

C/O 

Mg/Si 

Fe/Si 


Mass 

Mol 

Mass 

Mol 

Mass 

Mol 

Sun 

0.37 

0.50 

0.88 

1.02 

1.69 

0.86 

55 Cancri 

0.75 

1.00 

1.51 

1.74 

1.81 

0.91 

HD 967 

0.19 

0.25 

1.17 

1.35 

1.54 

0.78 

HD 1581 

0.32 

0.43 

0.97 

1.12 

2.39 

1.20 

HD 2071 

0.36 

0.48 

0.93 

1.07 

2.56 

1.29 

HD 4308 

0.29 

0.39 

1.14 

1.32 

1.77 

0.89 

HD 6348 

0.38 

0.50 

1.07 

1.23 

2.13 

1.07 

HD 9826 

0.47 

0.63 

1.02 

1.18 

1.72 

0.87 

HD 19034 

0.24 

0.32 

1.245 

1.45 

1.61 

0.81 

HD 37124 

0.31 

0.41 

1.46 

1.70 

1.34 

0.68 

HD 66428 

0.45 

0.60 

1.85 

2.14 

2.39 

1.20 

HD 95128 

0.41 

0.54 

0.86 

1.00 

1.73 

0.87 

HD 142022 

0.45 

0.60 

1.40 

1.62 

2.57 

1.29 

HD 147513 

0.12 

0.16 

0.66 

0.76 

2.18 

1.10 

HD 160691 

0.56 

0.74 

0.88 

1.02 

1.69 

0.85 

HD 213240 

0.34 

0.45 

1.27 

1.48 

1.99 

1.00 

HD 220367 

0.305 

0.41 

1.34 

1.55 

2.33 

1.17 

HD 221287 

0.24 

0.32 

0.75 

0.87 

2.81 

1.41 


of each population was computed along with the deviation 
from this value. The latter was used as the error bar. 


4.1. Mg/Si and Fe/Si ratios 

Figure ID provides mode values of Mg/Si and Figures 0 of 
Fe/Si in planets formed for each star composition versus the 
central ratio of the host star as given in Table Q] for a case 
with irradiation. For all three types of planets, and for the 
two irradiation cases, the element ratios between stars and 
planets correlate along a 1:1 relationship. This is an ex¬ 
pected result for icy and giant planets that mainly formed 
outside the ice line, a region where all refractory material 
has condensed. Inside the ice line, the elemental ratio in 
the disc is in some instances different from the stellar value, 
as shown in Figure [3j since the temperature is too high to 
condense all Mg- and Fe-bearing species. It is thus possible 
to form planets with non-stellar ratios, either because of 
their migration pathway or the location of formation, close 
to the star, which will result in different element ratios in 
the planets. However, the formation and migration of most 
rocky planets inside the ice line mainly occurs in regions 
where the temperature and the pressure during condensa¬ 
tion are such that these ratios will be stellar in the planets. 


4.2. C/O ratio in solids 

Figure 2] shows the mode value of C/O in solids (refrac¬ 
tory species and ices) for the three populations versus the 
stellar ratio. The C/O ratio in icy and giant planets ap¬ 
pears to be very weakly dependent on the C/O ratio in the 
star. In i- contras^o_other_studies_(see J _e i ^ iJ _Madhusudhan 
et al. 12012 h . an increase of the C/O ratio does not imply a 
carbon-rich interior in this study, which could be due to the 
assumption that the compos ition of the p lan etesim als does 
not evolvejas^^suegeste^^J^oriart^e^ayj^OlJ^Mori- 
arty et al. (|2014f ) also calculated successive condensation 
steps and their results show that successive condensations 
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Mg/Si planet (molar ratio) 

1 1.2 1.4 1.6 1.8 2 2.2 



in 

"So 


2.2 

2 

1.8 

1.6 

1.4 

1.2 

1 

0.8 


Fig. 1: Correlation of Mg/Si in solids accreted by the mod¬ 
elled planets and their respective host star, with irradiation. 
The line shows the 1:1 relationship. The results for simu¬ 
lations without irradiation are similar. 


Fe/Si planet (molar ratio) 


0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 



3 

<D 


Fig. 2: Same as Figure Q] for Fe/Si. 


favour the formation of C-rich material during later con¬ 
densations. The largest differences also come from the in¬ 
clusion of volatile species, which are not taken into account 
in these other studies, and which governs the amount of C 
incorporated into planetesimals. Another possibility is that 
there is not enough carbon to form such planets and that 
the C/O ratio of the host star needs to be higher for the 
typical pressures involved in protoplanetary discs. 

The difference between the types of planets can be ex¬ 
plained by their location of formation and migration path¬ 
way. Rocky planets are formed in a region where no ices 
are formed, thus the abundance of carbon is very low (see 
Fig. [3] - C is mainly used to form ices). Icy planets form 
beyond the ice lines, and the ices formed th ere control the 
C and O budget (see iMarboeuf et a.ljf2014bf) . This leads to 
a higher abundance of carbon and consequently to a higher 
C/O ratio. Giant planets fall between the two exemples. 
They are formed in the simulations mainly between the 
water and CO 2 ice lines, thus the C/O ratio in solids is 
lower than that of the icy planets, but higher than that of 
the rocky planets. For a model without irradiation, icy and 
giant planets have similar C/O ratios. In this case, ice lines 
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Fig. 3: Elemental mass ratios Fe/Si, Mg/Si, and C/0 
in the disc #1 of A13 (E 0 =95.8 g.cm -2 , a CO re=46 AU, 
7 =0.9), with solar composition and without irradiation. 
At a distance farther than 0.3 AU, the Fe/Si and Mg/Si 
ratios attain solar values (Fe/Sig=1.69; Mg/SiQ=0.88; 
C/O 0 =O.37). 


C/O planet (molar ratio) 

1 . 10’ 5 1 . 10' 4 1 . 10' 3 1 . 10' 2 1 . 10' 1 1 . 10 ° 



C/O planet (mass ratio) 


Fig. 4: Correlation of C/O in solids acrreted by the mod¬ 
elled planets and their respective host star, with irradiation. 
The results for simulations without irradiation are similar 
with values of roughly one order of magnitude lower. 


are closer to each other, which leads to minimal differences 
between icy and giant planets. Consequently, the C/O ra¬ 
tio in solids is mainly constrained by ices for non-terrestrial 
planets and is limited by the low amount of C in refractory 
species for terrestrial planets. 


Fig. 5: Time evolution of the nrid-plane temperature (solid 
lines) and pressure (dotted lines) of disc 1 (Eo=95.8 g.cm~ 2 , 
a C nro= 46 AU, y =0.9) in a case without irradiation. From 
iThiabaud et al.1 (I2015T ). 

and the C/OratiooftlmMiveloj^ofa^^lane^Jndeed^Thi- 
abaud et al. (120151 1 showed that two planets located at the 
same distance with the same mass could acquire a differ¬ 
ent C/O ratio depending on the time of accretion and their 
formation pathways. It is thus hardly possible to relate the 
C/O ratio in the e nvelope of a plane t to its host star. How¬ 
ever, as shown in IThiabaud et ahl (i2015h . the total C/O 
ratio (which includes ices and gas) is strongly dominated 
by ices. Consequently, the results of the total C/O ratio 
are similar to those presented in the previous section. 

4.4. Effect of T and P 

Mg/Si and Fe/Si ratios are always stellar beyond the water 
iceline since the temperature is low enough to condense 
all refractory material, which explains the absence of large 
error bars in Figures [T| and [2] for icy and giant planets. 
However, this is not the case at all locations inside the water 
iceline. The ratios at these heliocentric distances depend on 
what has condensed within the disc. Since thermodynamic 
properties are different among molecules, they condense at 
different temperatures and pressures, and thus the budget 
of Fe, Mg, and Si might differ depending on the temperature 
and pressure profiles of discs, resulting in non-stellar values 
for Mg/Si and Fe/Si ratios, as shown in Figure [3] Figure 
[5] shows the profiles of temperature and profile for disc 1 
of A13. Beyond 0.3 AU, the temperature is low enough (~ 
1000K) so that most of the Mg-, Si- and Fe-bearing species 
condense, which is not the case closer to the star. Hence 
ratios diverge from stellar values. 


4.3. C/O ratio in gas 

As discussed in IThiabaud et all (l2015f) , the dynamical pro¬ 
cesses involved in the C/O ratio of the envelopes of planets 
(migration, gas phase evolution, etc.) enables each type 
of planets to form gaseous envelopes with a range in their 
C/O mass ratios, ranging from 0 to 1.2 in the case of solar 
composition. Because of these processes, which do not de¬ 
pend on the chemical properties of the star, it is impossible 
to establish a relationship between the stellar composition 


4.5. In situ formation 

Migration is an important process that affects the chemical 
composition of a planet, but it can smoothen the compo¬ 
sitional differences of bodies formed at different distances 
to the host star. If planets are formed in situ, then the 
ratios (and especially the C/O ratio) will be representative 
of that ratio in a specific location of the disc. The effect 
of migration is not as strong on the Fe/Si and Mg/Si ra¬ 
tios, since the stellar value is easily obtained (see Fig.3) due 
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to the early condensation of these elements. However, it is 
stronger for the C/O ratio, because of the volatile character 
of most C-bearing species, as discussed in Section 4.2. 

Considering a giant planet formed at a final position of 5 
AU in disc 1 of A13. The water iceline in this disc is located 
at 1.516 AU, and all the other volatile molecules condense 
within 4 AU. In this case, the C/O ratio in the solids at 5 
AU will be at its maximum (all C and O has condensed, 
except for some losses during formation of volatile species), 
and thus the giant planets, which formed in situ, will have 
this ratio as well. If the giant planet migrated from 3 AU to 
5 AU, into a region where the C/O ratio is lower, a smaller 
amount of C present in the nebula has condensed compared 
to the amount of O because less C-bearing volatile species 
have condensed), then the C/O ratio at the end of the sim¬ 
ulation will be lower, resulting in a C/O ratio that would be 
representative of a mix between what was accreted between 
3AU and 5AU. Icy planets and rocky planets migrate less 
efficiently, resulting in higher (resp. lower) C/O ratios for 
icy planets (resp. rocky planets), where nearly all the C- 
and O-bearing species have condensed (resp. nearly none of 
the C-bearing species have condensed, but some O-bearing 
species have). 

Since the Mg/Si and Fe/Si ratios have attained a stellar 
value beyond 0.3 AU in this disc, the in situ formation 
of giant and icy planets will not change, except for hot 
Jupiters very close to their star. The change will be seen in 
rocky planets, however. If no migration occurs, the mode 
value of such ratios for rocky planets will not be stellar. 

5. Summary and conclusions 

The planet formation mod el of A13 was used alon g with 
the chemicalmodeljase^hi^arboeuf^tal. j2014aJ^_Mar- 
boeuf et al. (|2014bl l. and iThiabaud et all ( 20 1 4l . 1 201 oh to 
derive the composition of planets around 18 different stars 
with different Fe/Si, Mg/Si, and C/O ratios to deduce their 
elemental ratios. This study shows that the composition of 
planets that formed from condensed material from a neb¬ 
ula closely resembles the composition of the host star, as 
assumed in interior structure models, for solar masses and 
luminosities. 

The Fe/Si and Mg/Si ratios were found to be the same in 
planets and host star. The C/O ratio in solids was found to 
be very weakly dependent on the host star ratio for all pop¬ 
ulations, which contradicts previous studies on the carbon 
content of terrestrial planets, mainly because these studies 
did not consider volatile species. The total C/O ratio in 
planets is governed by ices, and thus the location of planet 
formation within the evolving planetary disc plays an im¬ 
portant role, as does possible migration of these bodies dur¬ 
ing and after accretion. The C/O acquired in the envelope 
of a planet can vary greatly between two populations due 
to physical processes, and it is thus impossible to link this 
ratio to the composition of the host star. This highlights 
how complex the origin of volatiles on rocky planets can 
be and emphasizes that it is necessary to study the origin 
of volatile components in our solar system by examining 
meteorites more close. 
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Thiabaud et al.: Elemental ratios in stars vs planets (RN) 


Table A.l: Elemental abundances (in mole) in stars selected for this sample study. The abundances of all elements in 
all stars are reported relative to an H abundance of 1 x 10 12 mole. When no values were found for an element, the solar 
elemental abu ndances were ta ken for t h e cal c ulatio ns. 

Reference s: 1. iLoddersI (120031) 2. Gilli et al. (|2006l) 3. Beirao et al.l (l2005li 4. Ecuvillon et al.| (12004 ) 5 . lEcuvillon et al.l 
(l2006fl 6 . iDelgado Mena et alJ (f201flh 7. iNeves et al.1 (l2009t) 8. ISousa et alJ (120081) 9. I Gonzalez fc Lawsl (120071) 10. Kang 

et al. (|201lh 11. iGaffau et al.l (I2005T ) 12. iTakeda fc Hondal (12005H 13. iMorel fc Micelal l|2004f) 


Star 

C 

O 

Mg 

Ni 

Si 

S 

Ti 

Fe 

Sun 

2.88 x 10 s 

5.75 x 10 8 

4.17 x 10 Y 

1.95 x 10 b 

4.07 x 10' 

1.82 x 10 Y 

1.00 x 10 b 

3.47 x 10 Y 

55 Cancri 

7.40 x 10 s 

7.40 x 10 8 

1.20 x 10 s 

3.60 x 10 6 

6.90 x 10 7 

2.10 x 10 7 

2.20 x 10 5 

6.30 x 10 7 

HD 967 

1.45 x 10 s 

5.75 x 10 8 

1.70 x 10 7 

3.98 x 10 5 

1.26 x 10 7 

- 

4.17 x 10 4 

9.77 x 10 6 

HD 1581 

2.29 x 10 s 

5.37 x 10 8 

2.88 x 10 7 

1.12 x 10 6 

2.57 x 10 7 

- 

7.41 x 10 4 

3.09 x 10 7 

HD 2071 

2.75 x 10 s 

5.75 x 10 8 

3.16 x 10 7 

1.38 x 10 6 

2.95 x 10 7 

- 

8.71 x 10 4 

3.80 x 10 7 

HD 4308 

2.63 x 10 s 

6.76 x 10 8 

3.16 x 10 7 

8.91 x 10 5 

2.40 x 10 7 

- 

7.59 x 10 4 

2.14 x 10 7 

HD 6348 

1.91 x 10 s 

3.80 x 10 8 

1.48 x 10 7 

4.57 x 10 5 

1.20 x 10 7 

- 

4.68 x 10 4 

1.29 x 10 7 

HD 9826 

4.79 x 10 s 

7.59 x 10 8 

6.38 x 10 7 

2.37 x 10 6 

5.42 x 10 7 

1.74 x 10 7 

1.35 x 10 5 

4.69 x 10 7 

HD 19034 

2.24 x 10 8 

6.92 x 10 8 

2.75 x 10 7 

6.31 x 10 5 

1.91 x 10 7 

- 

6.61 x 10 4 

1.55 x 10 7 

HD 37124 

2.46 x 10 8 

5.96 x 10 8 

3.24 x 10 7 

7.41 x 10 5 

1.91 x 10 7 

5.25 x 10 6 

7.00 x 10 4 

1.29 x 10 7 

HD 66428 

6.31 x 10 s 

1.05 x 10 9 

1.48 x 10 s 

3.72 x 10 6 

6.92 x 10 7 

- 

1.78 x 10 5 

8.32 x 10 7 

HD 95128 

3.84 x 10 s 

7.10 x 10 8 

4.27 x 10 7 

1.86 x 10 6 

4.27 x 10 7 

1.45 x 10 7 

1.14 x 10 5 

3.72 x 10 7 

HD 142022 

4.90 x 10 8 

8.20 x 10 8 

9.12 x 10 7 

2.82 x 10 6 

5.62 x 10 7 

- 

1.70 x 10 5 

7.24 x 10 7 

HD 147513 

2.40 x 10 8 

1.51 x 10 9 

3.47 x 10 7 

2.19 x 10 6 

4.57 x 10 7 

- 

1.62 x 10 5 

5.01 x 10 7 

HD 160691 

5.25 x 10 s 

7.08 x 10 8 

8.51 x 10 7 

4.37 x 10 6 

8.32 x 10 7 

2.57 x 10 7 

2.07 x 10 5 

7.08 x 10 7 

HD 213240 

5.40 x 10 s 

1.20 x 10 9 

6.50 x 10 7 

2.40 x 10 6 

4.40 x 10 7 

1.30 x 10 7 

1.30 x 10 5 

4.40 x 10 7 

HD 220367 

2.40 x 10 8 

5.89 x 10 s 

3.80 x 10 7 

1.05 x 10 6 

2.45 x 10 7 

- 

7.08 x 10 4 

2.88 x 10 7 

HD 221287 

3.09 x 10 8 

9.77 x 10 s 

3.16 x 10 7 

1.70 x 10 6 

3.63 x 10 7 

- 

1.12 x 10 5 

5.13 x 10 7 

Star 

N 

Na 

Al 

P 

Ca 

Cr 

Reference 

Sun 

7.94 x 10 Y 

2.34 x 10 b 

3.47 x 10 b 

3.47 x 10 b 

2.57 x 10 b 

5.25 x 10 b 

1 

55 Cancri 

1.80 x 10 8 

3.90 x 10 6 

8.70 x 10 6 

8.50 x 10 5 

2.80 x 10 6 

7.80 x 10 5 

2, 3, 

4,5 

HD 967 

- 

5.50 x 10 5 

1.23 x 10 6 

- 

7.94 x 10 5 

1.20 x 10 5 

00 

cd' 

HD 1581 

- 

1.38 x 10 6 

1.95 x 10 6 

- 

1.70 x 10 6 

3.31 x 10 5 

6, 7,8 

HD 2071 

- 

1.66 x 10 6 

2.45 x 10 6 

- 

1.95 x 10 6 

3.63 x 10 5 

6, 7,8 

HD 4308 

- 

1.17 x 10 6 

2.45 x 10 6 

- 

1.48 x 10 6 

2.14 x 10 5 

6, 7,8 

HD 6348 

- 

5.37 x 10 5 

1.17 x 10 6 

- 

7.94 x 10 5 

1.48 x 10 5 

6, 7,8 

HD 9826 

8.51 x 10 7 

3.94 x 10 6 

- 

- 

3.10 x 10 6 

- 

9, 11, 12 

HD 19034 

- 

9.12 x 10 5 

2.04 x 10 6 

- 

1.20 x 10 6 

1.78 x 10 5 

6, 7,8 

HD 37124 

- 

9.33 x 10 5 

2.29 x 10 6 

- 

1.15 x 10 6 

1.70 x 10 5 

9, 10, 11 

HD 66428 

- 

4.90 x 10 6 

6.31 x 10 6 

- 

3.72 x 10 6 

8.32 x 10 5 

6, 7,8 

HD 95128 

1.05 x 10 s 

2.40 x 10 6 

3.47 x 10 6 

- 

2.24 x 10 6 

4.90 x 10 5 

9, 10, 11, 12 

HD 142022 

- 

2.82 x 10 6 

5.13 x 10 6 

- 

3.39 x 10 6 

6.76 x 10 5 

6, 7,8 

HD 147513 

- 

1.95 x 10 6 

3.31 x 10 6 

- 

3.31 x 10 6 

4.90 x 10 5 

6, 7,8 

HD 160691 

- 

5.50 x 10 6 

7.59 x 10 6 

- 

4.47 x 10 6 

1.01 x 10 6 

2, 9, 

11, 13 

HD 213240 

1.40 x 10 7 

3.60 x 10 6 

4.70 x 10 6 

4.60 x 10 5 

2.50 x 10 6 

5.20 x 10 5 

2, 3, 

4,5 

HD 220367 

- 

1.29 x 10 6 

2.40 x 10 6 

- 

1.55 x 10 6 

2.88 x 10 5 

6, 7,8 

HD 221287 

- 

1.86 x 10 6 

2.40 x 10 6 

- 

2.63 x 10 6 

5.01 x 10 5 

6, 7,8 
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Table A.2: Characteristics of disc models from A13 assum¬ 
ing a gas-to-dust ratio of 100. 


Disc 

Mdisc (Mq) 

&core (AU) 

7 

1 

0.029 

46 

0.9 

2 

0.117 

127 

0.9 

3 

0.143 

198 

0.7 

4 

0.028 

126 

0.4 

5 

0.136 

80 

0.9 

6 

0.077 

153 

1.0 

7 

0.029 

33 

0.8 

8 

0.004 

20 

0.8 

9 

0.012 

26 

1.0 

10 

0.007 

26 

1.1 

11 

0.007 

38 

1.1 

12 

0.011 

14 

0.8 
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